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We employ first principles total energy and phonon calculations to address the structure and sta- 
bility of Al2Fe. This structure, which is reported as stable in the assessed Al-Fe phase diagram, is 
distinguished by an unusually low triclinic symmetry. The initial crystallographic structure determi- 
nation additionally featured an unusual hole large enough to accommodate an additional atom. Our 
calculations indicate the hole must be filled, but predict the triclinic structure is unstable relative 
to a simpler structure based on the prototype MoSi2 • This MoSi2 structure is interesting because it 
is predicted to be nonmagnetic, electrically insulating and high density, while the triclinic structure 
is magnetic, metallic and low density. We reconcile this seeming contradiction by demonstrating 
a high vibrational entropy that explains why the triclinic structure is stable at high temperatures. 
Finally, we note that AlsFe2 poses a similar problem of unexplained stability. 

PACS numbers: 61.66.-f,63.20.dk,64.30.-t,71.20.Lp 
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I. INTRODUCTION 

Aluminum-based intermetallic alloys with transition 
metals are of high interest for their complex crystalline 
and quasicrystalline structures, formed primarily with 
late transition metals, and their technologically useful 
compounds formed primarily with early transition met- 
als. Experimental phase diagram determination is diffi- 
cult because many phases often exist within small com- 
position ranges, many structures have unusually large 
unit cells and many are intrinsically disordered, exhibit- 
ing mixed or partial site occupancy. First principles cal- 
culations can help resolve some uncertainties in the phase 
diagrams, but are challenging themselves, for many of the 
same reasons. Intrinsic disorder requires studying alter- 
native realizations of specific site occupancy. Some of the 
nearby competing phases may have unknown or poorly 
known structures. The large unit cells pose computa- 
tional difficulties. Further complicating the study is the 
prevalence of magnetism among late transition metals. 

Although specific Al-Fe phases are not of direct com- 
mercial or technological interest (FIX THIS CLAIM!!!), 
precipitates of Al-Fe compounds can enhance the high 
temperature strength of pure Al. Al-Fe is also the pro- 
totype binary magnetic magnetic alloy based on a bcc 
structure [lj. Our own interest in the Al-Fe phase di- 
agram derives from its complex and disordered crystal 
structures, some of which are related to quasicrystals. 

The compound A^Fe is of special interest because of 
its unusual lowest-possible symmetry crystal structure, 
triclinic with space group # 1 (PI) as determined by 
Corby and Black [2J. The initial crystallographic refine- 
ment, based on anomalous dispersion experiments, pro- 
posed an 18-atom unit cell (Pearson symbol aP18) with 
an unusual "hole" that was sufficiently large to fit an 
entire Al or Fe atom. They also reported three sites of 



mixed occupancy, Alo.sFeo.s- 

Our preliminary first principles calculations of total 
energy [3] showed that filling the hole was energetically 
favorable, thus we predicted the correct Pearson type as 
aP19 3]. However, we found this structure to be unsta- 
ble with respect to competing phases of differing com- 
position, no matter how the hole was filled and how the 
partial occupancy was resolved. An alternate structure 
based on the MoSi2 prototype (Pearson symbol tI6) was 
predicted to be the true stable structure. This tI6 struc- 
ture can be considered as an Al-rich variant of the B2 
(Pearson cP2) structure of AlFe. It has never been ob- 
served experimentally, although it would be of high in- 
terest because it is predicted to be electrically insulat- 
ing with a narrow gap [J, lf| . Instead, multiple reexam- 
inations confirm the stability of a triclinic structure for 
A^Fe. A recent crystallographic refinement [6|], utiliz- 
ing conventional single crystal diffraction, proposes the 
space group is # 2 (PI) and fills the hole, confirming our 
predicted Pearson type aP19. 

Here we present a thorough study of the stability of 
A^Fe utilizing first principles total energy calculations 
of low temperature enthalpy supplemented by a phonon- 
based calculation of vibrational entropy yielding the high 
temperature Gibbs free energy. Our calculations predict 
that the tI6 structure loses stability to the aP19 structure 
at elevated temperatures. This occurs because aP19 has 
a much lower atomic density than tI6, resulting in high 
vibrational entropy. A similar stabilization effect due to 
vibrational entropy was observed in the 9/9' system of 
Al 2 Cu 0. 

The remainder of this introduction surveys the global 
Al-Fe phase diagram and presents our calculational meth- 
ods. We then present a thorough investigation of the en- 
ergetics of plausible aP18 and aP19 structures, including 
the effects of magnetic moment formation and antiferro- 



magnetism. Finally we present vibrational densities of 
states that display a large enhancement of low frequency 
phonons in the aP19 structure relative to tI6. 



A. Assessed Al-Fe Phase Diagram 

The Al-Fe phase diagram [8J contains at least six com- 
pounds as well as the two pure elements. Additionally 
there are at least three known metastable phases. TableU 
displays pertinent information including names, compo- 
sition ranges, Pearson types, space groups and assessed 
stability of all reported phases. Several of the phases re- 
port composition ranges associated with chemical substi- 
tution between Al and Fe and also partial site occupancy. 
The Al 3 Fe phase, more accurately described as Al 13 Fe4, 
is well known as a decagonal quasicrystal approximant. 
Structures of A^Fe and Al5Fe2 also feature pentagonal 
networks [9j. 
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TABLE I: Known phases of Al-Fe. Stability designation: 
S=stable to low T; S?=stable at high temperature down to 
unknown T<400C; M=metastable; HT=stable at high tem- 
perature only. See text section II Bl for AHf OT and A_E (units 
are meV/atom). 



B. Methods 

Our calculations follow methods outlined in a prior 
paper flOj. We utilize VASP [O, [13 to carry out first 
principles total energy calculations in the PW91 gener- 
alized gradient approximation. Comparisons with the 
LDA and PBE density functionals confirm the principal 
findings based on PW91. We relax all atomic positions 
and lattice parameters, and increase our k-point densi- 
ties until energies have conver ged. We adopt projector 
augmented wave potentials [13|, [l4 1 and maintain a fixed 
energy cutoff of 267.9 eV (the default for Fe). All cal- 



culations considered the possibility of spin polarization, 
and utilize the "Medium" precision setting which allows 
small wrap-around errors in Fourier Transforms. 

Given total energies for a variety of structures, we cal- 
culate the enthalpy of formation A_fff or which is the en- 
thalpy of the structure relative to a tie-line connecting 



the ground state configurations of the pure elements. The 
convex hull of these energies constitute the predicted low 
temperature stable structures. For structures that lie 
above the convex hull we calculate the instability energy 
AE as the enthalpy relative to the convex hull. When 
presented with mutiple structure possibilities, or mixed 
site occupancy, we examine all plausible structures and 
report the most energetically favorable. 



II. RESULTS 
A. Global phase diagram 

As illustrated in Fig. Q] our calculated enthalpies agree 
in almost every respect with the assessed phase diagram. 
All the phases that are known to be stable at low temper- 
ature (heavy circles) indeed reach the convex hull. Those 
known to be metastable (diamonds) or stable only at 
high temperatures (light circles) lie slightly above the 
convex hull. Of those whose low temperature stability 
is uncertain (triangles), A^Fe lies on the convex hull, 
while Al 5 Fe 2 and Al 2 Fe (aP18 and aP19) both lie slightly 
above. The only serious descrepancy between the experi- 
mental phase diagram and our calculation is the presence 
of A^Fe in the tI6 (MoSi2) structure on the convex hull. 

Magnetism was found to be favorable in all the struc- 
tures reported in Table U containing 33% or more Fe. 
Elemental Fe exhibits fcrromagnetism in its low temper- 
ature cI2 structure and modulated antiferromagnetism 
in its high temperature cF4 structure. A^Fe also ex- 
hibits long-wavelength antiferromagnetism in its aP18 
and aP19 structures. 

Our calculated total energies correctly predict the 
AlgFe 2 and Al 6 Fe structures to be metastable in the Al- 
Fe alloy system, while the same structures are correctly 
predicted to be stable in the Al-Co and Al-Mn alloy sys- 
tems, respectively [3j, |l5[. Similarly, we correctly predict 
Al-Fe to be unstable when placed in the AlnMn4.aP15 
and AI5C02.I1P28 structures. Conversely, we correctly 
predict AlFe3 stable in Al-Fe but unstable in both Al- 
Co and Al-Mn. This sensitivity to the small differences 
in interatomic bonding between Fe and its neighbors in 
the periodic table, Co and Mn, gives us confidence in the 
validity of our first-principles total energies. 



B. Al 2 Fe.aP18/aP19 

The Corby and Black [2J structure with Pearson type 
aP18 possesses an unusually large hole. In addition three 
sites, labeled Ml, M2 and M3, exhibit mixed Al . 5 Fe .5 
occupancy. We tested all eight arrangements of Al and 
Fe among these three sites (labeled AAA, AAF, . . . , 
FFF) within a single unit cell and found that none of 
them resulted in a stable structure. Additionally, all 
suffered rather large maximum atomic displacements of 
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FIG. 1: Enthalpies of formation of Al-Fe compounds. Plotting symbols: heavy circle=known stable; light circle=known high 
temperature; diamond=known metastable; triangle=unknown stability; square=unknown phase. Plotting colors: black=on 
convex hull; red=above convex hull. Notation A and F indicate Al or Fe on Wyckoff la positions Ml, M2 and M3 (aP18) or 
Wyckoff 2i position Fe4 (aP19). 



0.3-0.4 A. All structures except AAA favored weak ferro- 
magnetism. Our optimal stoichiometric structure, AFA, 
exhibited magnetic moments averaging 1 fj,g/Fe atom, 
resulting in an energy drop of 26 meV/atom. 

Previously Q, we found that filling the Corby-Black 
hole with an Fe atom is energetically favorable, sug- 
gesting that aP19 is the correct Pearson type for this 
compound. The nominal composition A^Fe cannot be 
achieved in a single unit cell. Additionally, in the new 
structure refinenement [g] there is a single site (Fe4, 
Wyckoff position 2i) with mixed Alo.295Feo.705 occu- 
pancy. For this reason we examine both Al and Fe on the 
Fe4 site within a single unit cell (labeled AA, AF=FA, 
FF), and we also examine a supercell doubled along the 
a-axis with one of the four Fe4 sites replaced by Al (la- 
beled FFFA). 

In our single cell aP19/FF ferromagnetic structure we 
find moments of: 1.1/ib on Fel; 0.9 on Fe2 sites; 0.8 
on Fe3; 1.8 on Fe4. In the supercell we find long-period 
antiferromagnetism is favored over ferromagnetism by 3 
meV/atom, consistent with the findings of incommen- 
surate antiferromagnetism [la ]. Magnetic moments in 
our aP19/FFFA antiferromagnetic structure were quali- 
tatively similar in magnitude to our aP19/FF ferromag- 
netic structure except in the vicinity of Al substitutions 
on the Fe4 site. Attempts to find disordered magnetic 
ground states by starting from randomly selected initial 



Fe moments of ±2/ib did not succeed in lowering the to- 
tal energy, suggesting a possible absence of spin glass or- 
der [171 ]. Likewise, runs utilizing noncollinear magnetism 
failed to reduce the energy. 

We note a slope in the enthalpies of aP19 with respect 
to XFe so that AE drops to as low as 9 meV/atom in the 
Al-rich limit of aP19/AA (not shown) at x Fc = 0.2632, 
while it rises to 33.3 meV/atom in the Fe-rich limit of 
aP19/FF (reference number 12 in Fig. [J) at a; Fc =0.3684. 
The optimal structure we found at the precise xp c — 
0.3333 stoichiometry of ALjFe is a six-fold supercell of 
aP19 (not shown) containing 114 atoms and remaining 
unstable by 29 meV/atom. The experimentally reported 
composition is a;F e =0.3374. 



III. VIBRATIONAL AND CONFIGURATIONAL 
FREE ENERGIES 

Despite extensive effort, we were unable to find an 
enthalpy-minimizing triclinic structure for A^Fe. All 
structures considered remained above the convex hull as 
illustrated in Fig[TJ Given that the assessed phase dia- 
gram does not assert low temperature stability, we now 
consider the possibility that this structures is stabilized 
at high temperature by entropy. An entropy s per atom 
results in a free energy reduction of Ts at temperature 





FIG. 2: (color online) tI6 and aP19/FFFA vibrational (left) and electronic (right) densities of states. For convenience the 
electronic DOS of aP19 is shown without spin polarization. Vertical dashed lines indicate Fermi energy. 



T. The two primary sources of entropy are continuous 
vibrational and discrete configurational degrees of free- 
dom. 

Even with the hole filled, the atomic volume of 
aP19/FFFA is 13.48 A 3 /atom while the atomic volume 
of tI6 is only 12.7 A 3 /atom. The higher atomic volume 
of aP19 implies weaker bonding, consistent with lower 
enthalpy of formation. Weaker bonding suggests an en- 
hanced density of low frequency phonons that can in- 
crease vibrational entropy and correspondingly lower the 
vibrational free energy. Because Al2Fe.aP19/FFFA is at 
very nearly the identical composition to A^Fe.tK, we 
need only compare the vibrational free energies of these 
two structures (in principle we should include a small 
admixture of AlFe.cI2, but with very low weighting). 

The vibrational free energy can be obtained within the 
harmonic approximation, given a vibrational density of 
states g{ui), by means of the quantum mechanical theory 
of phonons as expressed by 



/, 



vib 



(T) = k B T J g(u)\n[2smh(huj/2k B T)}du 



(1) 



To calculate g(ui) we employed the force constant method 
for phonon calculations [18[ , in which we take a supercell 
of minimum edge length 8 A and evaluate the second 
derivatives of total energy <9 2 £//dRi<9Rj for all pairs of 
atoms i and j. The calculations utilized density func- 
tional perturbation theory as implemented in VASP 5.2. 
We then construct the dynamical matrix [19] D(k) on a 
dense mesh of fc-points and evaluate the phonon frequen- 
cies. For our phonon calculations we revise our standard 
VASP calculations to use "Accurate" precision so as to 
avoid wrap-around errors, and we set EDIFF=5xlO -10 . 
The resulting calculated vibrational densities of states 
are presented in Fig. [2] The excess density of low fre- 
quency phonons is clearly evident and yields a reduction 
in free energy of 38 meV/atom by the time we reach 
T=1000K. 



In addition to the continuous vibrational degrees of 
freedom, some entropy is contained within the discrete 
configurational fluctuations associated with the partially 
occupied state. We estimate this assuming independent 
occupancies of every mixed or partially occupied site. For 
the case of AFjFe, the mixed Al/Fe occupancy on the Fe4 
site yields 
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(2) 



where the factor of 2/19 arises because there are 2 Fe4 
sites/cell and the occupancy factors are Jai =0.295 and 
f Fe = 1 - fAl=0.705. Setting T=1000K we find that 
mixed occupancy can account for up to 5 meV/atom of 
stability. This is an upper bound on the actual value, 
because correlations in site occupancy will slightly reduce 
the entropy. 

Combining the vibrational free energy and the addi- 
tional discrete configurational entropy results in a signif- 
icant lowering of the Gibbs free energy of Al2Fe.aP19. 
We find this free energy reaches the convex hull around 
T=500K. The competing tI6 phase loses stability by 
T=550K. 



IV. AL5FE2 

The other phase of problematic stability is AlsFe2 (see 
Fig. [3]) . As noted previously this structure features pen- 
tagonal networks as well as "aluminum channels" . Alu- 
minum atoms can diffuse along the channels in a liquid- 
like fashion, and the density of the channels (i.e. the 
number of aluminums per c-repeat per channel) can vary 
over a range around one Al/repeat/channel. We find the 
energetically optimal structure occurs with four Al atoms 
in three repeats for a channel density of 4/3. The addi- 
tional Al atoms raise the Fermi energy to the vicinity 
of a strong pseudogap in the electronic density of states, 



which can account for its energetic optimality. However, 
the energy remains high relative to the convex hull by 
about 6 meV/atom. Preliminary phonon calculations fail 
to explain the observed high temperature stability of the 
phase. 

We believe that full anharmonic vibrational free energy 
calculations are needed, owing to the liquid-like motion 
of the channel Al atoms. This behavior is evident in 
Fig. [3] where we present the occupation densities p a { r ) 
using a color scheme where gray represents low density, 
while red and green represent high density of iron and 
aluminum, respectively. The data was collected using 
VASP molecular dynamics for a duration of 15ps. 



clinic structure aP19 when a simpler and energetically 
preferable tetragonal tI6 structure is available. Since 
aP19 is stabilized at high temperature by its high vibra- 
tional and discrete configurational entropy, it will natu- 
rally form first in samples cooled from the melt. Other 
high temperature phases at nearby compositions such as 
AlsFe2 and AlgFes compete with A^Fe, likely further 
frustrating attempts to achieve the tI6 structure through 
off-stoichiometric growth. Still, density functional theory 
predicts the true low temperature structure is tI6, and we 
encourage experimentalists to search for it. The payoff 
would be obtaining an insulating compound consisting of 
metals Al and Fc. 



CONCLUSION 
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FIG. 3: (color online) Structure of AlsFe2 showing aluminum channels. The 6-axis is horizontal, and the c-axis is vertical, (left) 
ideal period 3 superstructure, (right) Ab-initio MD at T=1300K. 



